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Abstract 



We introduce a new class of spatial-temporal point processes based on 
Voronoi tessellations. At each step of such a process, a point is chosen at 
random according to a distribution determined by the associated Voronoi 
cells. The point is then removed, and a new random point is added to 
the configuration. The dynamics are simple and intuitive and could be 
applied to modeling natural phenomena. We prove ergodicity of these 
^ ■ processes under wide conditions. 
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1 Introduction 
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A spatial point process is a stochastic model for the location of events in a 
space, and as such it is a random element which takes discrete sets of points as 
its values. Trees in a forest, schools in a city, or capillaries on the surface of a 
bodily organ are examples of situations that can be modelled by a point process 
over a two-dimensional manifold, impurities in metals or positions of submarines 
require models denned over three-dimensional manifolds. The one-dimensional 
point process is often used as a model for events in time, but may also be applied 
to such problems as the physical distribution of files on a hard-disk where the 
underlying space can be treated as an interval in R. 

The Voronoi tessellation is a useful tool for the analysis of spatial point 
processes. This can be defined for any locally finite set of points A in a metric 
space (M, d) . With each x G A we associate its Voronoi cell 

C*{M) :={yeM: d(y } x) = min [d(y, x') : x' 6 A]} . (1) 
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[We will generally suppress the reference to M in the notation for the cell and 
just write where the underlying set is clear from the context.] The tessella- 
tion produced by X is 2?x = {C* : x G X}, and the elements of X are called 
the generators of 3?x- Statistics drawn from analysis of this structure provide 
an intrinsic description of the distribution pattern of the set of generators and 
have accordingly been investigated extensively (see e.g. pQ). The Voronoi tes- 
sellation is also a natural object of interest whenever concepts such as catchment 
area or zone of influence are appropriate to the situation being modelled, as in 
most of the examples mentioned above. When this is the case it is also very 
interesting to consider point processes evolving in time under laws which are 
determined by characteristics of the Voronoi tessellation. A simple example is 
a process in which the lifetime of a point is a random variable with distribution 
determined by its Voronoi cell. Examples of this kind of spatial-temporal point 
process are the Hotelling processes [13], which model the geographical distri- 
bution of businesses competeing for market share and the Adjustment models 
for territorial animal behaviour ^3]. As far as we know there has been no gen- 
eral treatment of such models although they promise to produce a rich class of 
stochastic objects typifying various kinds of spatial clustering. 

In this paper we consider some classes of such models viewed as discrete- 
time Markov processes {X n } n>0 taking values in a fixed finite-dimensional con- 
figuration space. Let (M, d) be a metric space and N a fixed positive integer 
representing the initial number of points in M. The associated configuration 
space will be either i£T = (J fc=1 M k or SC = M N , depending on whether the 
process is of a thinning or thinning-replacement type. Any finite set of points 
X generating a Voronoi tessellation can of course be represented in many dif- 
ferent ways as a point x e 5£ since the order of the points is immaterial to the 
tessellation, however it is convenient to retain an ordering of the points in the 
configuration space. In the thinning-replacement process the total number of 
points after each replacement remains constant, although the thinning and re- 
placement components are independent stochastic processes. Replacement this 
will always be determined by a probability measure /i on M, equivalent to the 
volume measure, whereas the thinning follows a probability rule on the Voronoi 
tessellation. We tacitly assume that all subsets of M or iT appearing in this 
paper are measurable. 

Let a selection function S : ^ — > M + be given, defined on the space ^ of 
all possible Voronoi cells in M. If x = (x%, . . . , x^) is the current configuration, 
then in the next step exactly one co-ordinate point Xj, 1 < J < k, is cho- 
sen at random and either removed (in a thinning process) or reassigned (in a 
thinning-replacement process) to a random position. The probability of choos- 
ing coordinate Xj is proportional to the value of the selection function S(C^ ) 
on the respective Voronoi cell C*., i.e. 



3(Cl) 

p r fj — j\ x ) •— v x i' (2) 

There are many quantifiable properties of a Voronoi cell upon which the func- 
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tion S can be based. Some of these such as volume, perimeter or surface area 
or suitable high dimensional generalisation, number of edges, faces, minimal or 
maximal internal angles etc., are properties of any closed simplicial complex in 
the same underlying space. Others, such as volume of the associated Voronoi 

flower (Voronoi flower of C%. is the closure of the set jz G M : Cff^ ^ Cf. j) 
or the distance to nearest neighbour, are specific to the Voronoi structure but 
can be computed from the properties of the individual cell. A further generalisa- 
tion could allow S to be a function of Xj and of x, for example the weighted sum 
of cell functions in the original sense with weights determined by the number of 
edges in the smallest nearest-neighbour arc to Xj. In this sense, the degree of a 
selection function S could be defined as the depth of the Voronoi nearest neigh- 
bour relation required to fix its value. We will consider only the simplest kinds 
of first-degree selection functions. More specifically, we consider the following 
two classes of such processes. 

(A) The "volume-based" , or t>-process. We require enough additional struc- 
ture for (M, d) so that cell-volume can be defined. Let M be a compact 
piecewise-smooth manifold equipped with a measure A, which is equivalent to 
the volume measure on M with a density bounded away from zero and from 
infinity. We assume that the value of S(C£.) is determined by the value of 
Air;''): 

s{c: 3 ) = s a \i( ':;:)). 

Without ambiguity we drop the subscript in S v and consider S : 1R + — > M + to 
be a function of A(C* ). If S is increasing, then points with Voronoi cells of 
large volume are more likely to be chosen to be culled or moved, and so the 
selection pressure favours points with small cells, that is, points restricted by 
close neighbours. A decreasing S favours points with large cells, that is, isolated 
points or points whose near neighbours fall within a limited arc. Functions of 
the form S(v) = v a , a G K. produce scale-independent models. 

(B) The neighbour-based, or n-process. The Voronoi tessellation determines 
for each generator a set of its Voronoi-nearest neighbours. Formally, this set is 
defined as 

[Xj] x : = |xi : x = Oi, ...,x k ),i^j, card(C* n C* ) > 1 j . 

We now assume that S(C£.) = S n (ca.id([xj] x )), for some S n : {1, . . . , N} R, 
where again we drop the subscript where the context is clear. The n-process 
requires less structure on (M, d), but the selection function S determines which, 
if any, types of cells are favoured by the evolution. 

Figures^ to 0] below illustrate some of the behaviours that were observed in 
simulations. 

Fig. ^depicts side-by-side realisations of three different f-processes on a circle 
with S(v) = v a , each having the same total number of points (N = 128), but 
different a values. The base of each rectangle represents the circle M opened 
out into a line segment by a cut, whereas the y-axis represents the time. The 
well defined clustering observed in (c) was found to occur when a > 1 and iV 
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Figure 1: Evolution of w-processes on circle. Number of points, N = 128, number of steps, 
T = 4096, S{v) = v a , with (a) a = -1.0, (b) a = 0.5, (c) a = 1.5. 

is sufficiently large. The phase change was observable even for values of a close 
to 1, as seen in Fig.El where the time scale has been compressed by a factor of 
5 from Fig.^ to bring it out more clearly. 

Fig. El shows the results of running four instances of the f-process on the 
unit square with S(v ) = v a and values of a ranging from -3.0 to 1.4. The phase 
change at a = 1 was observed just as in the one-dimensional case. For values 
of a < 1 a smooth gradation of degrees of clustering was observed without 
apparent interference from edge effects. 

Two instances of the n-process are shown in Fig.0] The n-processes pro- 
duced a very rich collection of different clustering behaviours, but in general it 
would seem that selections which favour cells with greater than average num- 
bers of neighbours (average is six) lead to less clustered configurations, and vice 
versa. 
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Figure 2: Phase change in the w-process on circle. N = 128, T = 20480, S(v) = v a , (a) 
a = 0.95, (b) a = 1.0, (c) a = 1.05. 
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Figure 3: Realisations of the w-proccss with S(v) = v a for -3.0 < a < 1.4. Left-Right,Top- 
Bottom: a =-3.0, 0.2, 1.0, 1.4. In each case there are N = 2000 points with the number of 
cycles T = 24000. 
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Figure 4: Realisations of the n-process with N = 2000, T = 16000. The pattern on the 
left was produced by the selection function S(n) = 0.1 + \n — 6| 2 which "favours" cells with 
six or close to six sides, that on the right by S(n) = l,n/ 5;S'(5) = 5000 which strongly 
"disfavours" cells with five sides. Note the apparent edge-effects in the former case. 
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We will discuss the simulation results in more detail, and consider the prob- 
lem of inferring the selection function from the statistics of a given point pattern 
which has arisen from a v- or n-process elsewhere. In particular, we employ 
Thiel's redundancy measure pQ and Baddeley and Van Lieshout's J-function jlj 
as a first step in the characterisation of these patterns. 

Simulations appear to demonstrate ergodic-type behaviour for t>-process 
with S(v) = v a when a < 1 the behaviour has a very different character which 
indicates that ergodicity may not take place. In the present paper we solve the 
problem of ergodicity raised by these observations. 

In Section 2 we prove, under rather broad conditions, the ergodicity for 
both the n- and v-processes. In Section 3 we present a crude explanatory model 
which explains the phase change that is observed at a = 1 in the w-process with 
S(v) = v a . 

Remark 1. Our process with replacement resembles a Gibbs sampler in that at 
each step only one coordinate of the state is resampled. The main difference 
is that in a Gibbs sampler the choice of the coordinate to be resampled is 
deterministic, but the new sample is from a distribution conditioned on the 
other current points. In our case, the choice of point is from a distribution 
conditioned by the set of current points, and the new sample is i.i.d. This latter 
choice could easily be modified so that the new positions of the resampled 
coordinates are governed by a general set of conditional distributions, as in the 
usual application of the Gibbs sampler. 

2 Main results 

In this section we consider only processes with replacement. Assume M is 
a compact piecewise-C 2 manifold in R m , m > 1, with or without boundary, 
endowed with the geodesic metric d. Let fi be a probabiloty and A a measure 
on M, both equivalent to the volume measure v on the manifold. 

We define a discrete-time Markov n- or f-process with state space SC = M N 
as follows. Let x G 3C , A a Borel subset of JT, z G M, A C M, 

x \ji z ] := • • • > x j-li z i x j+lj • • • > x N)i 

Aj[x] := {z G M : x[j, z] G A} and 
A[j, A] := {x[j, z] G : x G A, z G A} . 

We also write, with a slight abuse of notation, x G x when x = (x\, . . . , xn) 
and x = Xj for some j, 1 < j < N. 

The transition probability function is given by 

N 

P(x,A) := Pr(X n+1 G A\X n = x) = ^ ^{A 3 [x])Yi{ J = j\x), (3) 

3=1 

where Pr(J = j\x) is the selection probability function @. 
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Theorem 1. (i) The n-process with replacement with a selection function S : 
{1, . . . iV} — > IR + which takes only positive values is Harris ergodic. 

(ii) The v-process with replacement with a selection function S : (0, A(M)] — > 
M. + , such that both S and 1/S are locally bounded, is Harris ergodic. 

Remark 2. It's not difficult to see that the assertion of part (i) holds under more 
general conditions as well. 

Proof, (i) Denote by fi® N the iV-fold product measure on S*, and by P^\ the 
/c-step transition function generated by P. We start with the following lemma. 

Lemma 1. The n-process with replacement, with M and fi as specified above, 
satisfies the following Doblin condition: there exist e < 1,7 > such that if 
H® N {A) > e for a measurable A C 3£ ' , then P^ N \x, A) > 7 for any x G 3£ . 

Proof. Let B(z, 5) = {x e M : d(x, z) < 5} be the 5-ball around the point z 6 

N 

M, and let H — \\ B(zi,5i) C SC ', for z { e M,5i > 0. It suffices to find a 

i=l 

constant k > 0, independent of {z i: 5i : 1 < i < N}, such that for any x G 3£ 

P( n \x,H) > k ^ N (H). (4) 

Indeed, since any open subset of can be approximated in //^-measure ar- 
bitrarily well by disjoint unions of basic sets H, then by the regularity of the 
measures on E£ the relation Q implies that 

p( N \x,A)>k ^ N (A) (5) 

for any x G 3£ and A C 3C '. 

First, note that for any x = (x\, X2, ■ ■ ■ , xpj) G 3£ and A C M by @ and 

P(x,{x}\j,A]) >MA), (6) 

where 

k= A min ^^, >0. (7) 

Therefore, for any y G M and i ^ j we also have 

P{x[i,y),{x[i,y)}[j,A))>kfx{A). 

Taking into account that the Markov chain can proceed from x to the set 
{x} [i, B(zi,5i)][j, B(zj,5j)] by first changing the value of the i th component 
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and then that of the j th one or vice versa, we obtain from © that 

P^ (x, {x} [i, BM)U BizjM) 

P(x,{x} \j,dy j ))P(x\j,y j ),{x} \j,y,][i,dyi]) 



B( Zi ,8i)B(zj,8j) 



+ J J P(x,{x}[i,dy i ])P(x[i,y i \,{x}[i,y i \[j,dy j }) 



B(zj,8j) B(zi,8i) 
,2 



>2k 2 j nidyi) J ti(dy j ) = 2k 2 ti(B(z i ,5 i )) f i(B(z j ,5 j )). 
Continuing in this way, we will clearly get 

N N 

(x, B( Zj , 5,)) > N\k N ft ^B( Zj , 5,)), 

3=1 3=1 

which proves (j3J) and the lemma. 

□ 

The relation (JSJ clearly implies the irreducibilty of the n-process. Indeed, 
© also implies that for any x 6 , A C JT, 

p( N+1 \x,A) > k/j l ® N (A), 

and hence the process is aperiodic. Now the desired assertion is immediate from 
a theorem in (D . p. 391. Part (i) of the theorem is proved. 

(ii) Ergodicity of the f-process requires more analysis. We wish to produce 
an argument similar to that of Lemma 1. The measure \x determines probabil- 
ities for points entering the configuration (the replacement points) and it will 
be sufficient for our arguments to note that for any fixed S, 

Us ■= inf {fi(B(x, 5)) : x G M} > 

as a consequence of the smoothness and compactness of M. Similarly we can 
define, 

\ 5 := M{X(B(x,5)) : x G M } > 0. 

For the v -process we do not have a simple analogy to the lower bound derived 
from (J7J). Here what is required is a positive lower bound for X(B(x, 8) fl C*), 
which does not exist if x is unrestricted. For example, if x G dM, then a single 
other point in the configuration x can cause the Voronoi cell around x to be 
arbitrarily small in volume and, consequently, the value of the selection function 
may become arbirarily small on this cell. So first we need to define those points 
which are not too close to the boundary. 
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Let M s := {x G M : d(x, y) > 5 for any y G dM} for a 5 > 0; if M is without 
a boundary, then M5 := M. The smoothness of M guarantees that fi(Ms) > 
and X(Mg) > for 5 small enough. Take any x G M<5 and a; G and define 

, J card({y Ga; :0< d(x,y) < 25}) for rr G a;; 

\ iV otherwise. 

So if a? is such that n(x, £E,5) < 1, then B(x,5) intersects with at most two of 
the Voronoi cells C% . This means that if x G Mx is fixed, then for all 5 > 
sufficiently small, and any x G 3£ such that n(x,x,S) < 1, we have, via the 
smoothness of M, 

\(B(x,5)nC:) > X 5/3 . (8) 

Define for 5 > 0, 

@s '■= {( x i: ■ ■ ■ i x n) £ '■ all Xi G M,5 and Xj) > 25, i 7^ j} . 

Clearly we can choose 5 small enough so that fi® N \^s) > 0. 
Now define a new measure on i?f as follows: 

0(A) :=^(in^). 

Lemma 2. (i) There are 5 > and e > snc/i i/iat /or every x G , 

P W (^,)>e. (9) 

(ii) T/iere is a 7 > such that for 5 as above and any x G 3>s, and any Borel 
subset AC f 

P( N \ X ,A) > 70(A). (10) 

Proof, (i) Let 5o > be small enough that Ms has postive /x (and hence positive 
A) measure. By appealing to the smoothness structure of M or otherwise we 
can find N + 1 disjoint open subsets of Ms . Denote these by L\, L 2 , . . . , L N+1 
and fix them for the rest of the proof. Let 

5i := sup {5 : there are N disjoint balls B(zk, 35) C Li}, % — 1, . . . , N + 1. 

Clearly, 

5 := min <5j > 0. 

0<i<JV+l 

Now for any x = (x±, ...,%) G at least one of the Lj will be disjoint from 
{x±, . . . , Xn}- Therefore for such an % we can find \ . . . , z$ G Lj such that 

B(zf,35) n 5(4°, 35) = for j ^ fc, 

x k i B{zf,35) for all k,je{l,...,N}. 

Obviously, Uf=iB(zf\5) C 0*. 

As in Lemma 1, we consider a possible sequence of N steps of the Markov 
process in each of which exactly one of the co-ordinate points x m of x moves 
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into a unique B(zfi, 8),j — l,...,N. At each step, the probability of the chosen 
point landing in a B(z^) , 5) is greater than or equal to p,s, so we need only be 



concerned with the probability of all N original co-ordinate points moving in 
N consecutive steps of the process. This will have a lower bound independent 
of the initial state x if we can find a bound k > with the following property: 
whenever we are at a stage at which some, but not all, of the x m have moved 
and are in corresponding sets B(z^), 5), there is a j such that Xj has not moved 
yet, and 

Pr(J = j\x') > k, where x' is the current state of the process. 

Let us assume that in the first k > 1 steps k different original points 
, . . . , Xi k have moved. We denote their new positions by x f tl , . . . , x\ k . Write 
v :— {ii, . . . , ik}, £ := {1, . . . , N} \u = {j : Xj has not moved} and 

N 

■.= s{\(cZ)) = + E = : s * + & 

i=i j€" je£ 

Thenp := S|/S s = Pr(J G f| x') is the probability that at the next step one of 
the remaining points will move. 
If p > 1/2, then for some j G £, 

1 1 

Pr(J=7'|0> ttt > , (H) 

V Jl ' ~ 2card(0 ~~ 2N K J 

otherwise 

Si < S£. (12) 
We only have to consider the latter case. For an x'j, j G u, we have 

4 G 5(^,5) 
for some j', and so not only do we have 

d(4,<J>2<5, 

when j ^ m G u, but also for m G £, since 35) was free of any rc^ from 

the original state x. This implies that for the set W := Ujeg B(xj, 5) one has 
C5' fl W — for i G v, and therefore 



Hence, 



max A(Cr.) > 77— > 7--- > — , 

jet 3 ~ card(£) ~ card(£) ~ N 
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while for all j £ u, 

A(C^) > \s- 
Due to the local boundedness of S and S^ 1 , 

0<b:= inf S{v) < sup =: B < oo. 

A«/iV<u<A(Af) A 4 /At<d<A(M) 

Therefore by and ((EJ) 

b b 

maxPiU = j\x') > — > — — =: « > (13) 

i«* V ^ " S| + card(z/) J B " V ; 

Clearly k < so from (jTTj) and (fT5j) we have, by a similar argument to that 
in the proof Lemma 1, that 

pW{x,@ s )>N\k n - 1 ^, (14) 

which proves 0. 

(ii) Let Ac f and assume that 4>(A) > 0. As in the proof of Lemma 1, it 
will suffice to restrict attention to basic open sets 

N 



H := Y[ B{ Wi , n) C 3> s , Wi G M s 



8=1 



and as in the proof for (i), we specify a possible sequence of N moves in the 
Markov process and aim to produce a bound analogous to (fT^j) for P^ N \x, H). 
The term /if in (JUjl is easily seen to be replaceable in that case by 



N 



Y[ fi (B(w i ,r i )=<f>(H), 



i=l 

since it represents the product of the probabilities that the points move into the 
target sets; the value k from (i) is no longer valid however, because, retaining 
the notation above, we don't necessarily have that d(xj, x' m ) > 25 when j £ £ 
and m £ v. Thus, if we can find a new constant k' > such that at every 
intermediate state x' 

maxPr(J = j\x') > k', (15) 

we will have 

P( n \x, H) > NW N - l (j>(H), (16) 

and ()1U|) will have been proven. 

As before, let x' be the state of the process after k steps, and let j £ £. 
Since x £ for any other z £ £, we have d(xj,Xi) > 25. On the other hand, 
for all Z,m £ v we also have ^(x^x^) > 25 as these points will be part of the 
configuration that we will get after N steps, which has to belong to H C @$. 
Let m £ f be such that 

(i(x^,Xj) = mm.d{x^Xj). 
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If d(x' m ,Xj) < 5/2 then for all other i G v, d(x' i ,Xj) > 35/2, and so the ball 
B(xj,35/A) is contained entirely within C*. U C*, , and hence for sufficiently 
small 5 > 0, by ©, 

A(C*.) > Atf/ 8 . 

Alternately, d(x' i ,Xj) > 5/2 for all i G z/, and so B(xj,6/A) C (7*'. In any case, 
A(C^') > A5/8, and hence, 

A5/8 



Pr(J = j|aj')> <, 



We also have, 



Ss < card(£) sup S'(f) + card(z/) sup S^v) 

«e[Aj /8 ,A(Af)] oe[Aj,A(Af)] 

< N sup ^(i;) =: Q < 00. 
ve[\ s/8 ,\(M)] 

So we can take k' := Q^ 1 ini v€ [x s/s ,x( M )} ^( v ) i n (ffSJ- Lemma 2 is proved. □ 
Following § 5.4.3 of [J2|, aperiodicity for the v-process is defined as follows. 

Let 

E» 6 := {n > 1 : there is a 7 n > s.t. P {n \x, A) > -f n (f)(A) 

for any x G @ s , A C JT} . 

The u-process is called aperiodic if g.c.&.{E@ s ) = 1. 
Lemma 3. JTie v-process is aperiodic. 

Proof. It is enough to find a 77V+1 > such that for any x G ^5, and any Borel 
subset AC J%~ 

p( N+1 \x,A)> lN+1 <P(A). (17) 

From Lemma 2(ii) 
p( N+1 \x,A) = [ P(x,dy)p( N \y,A) 

> [ P(x,dy)p( N \y,A) > P(x,@ s ) 1( f)(A). 

JS>5 

Without loss of generality we can choose 5 so small that for any N points 
Zi G Ms, 1 < i < N, such that the balls B(zi,8) are disjoint, we can find a 
zn + i G Ms with 

B(z N+1 ,S)n |J = 0. 

l<i<iV 

For such a 5 it is clear from (jSJ) and the definition of that ^5) > //«s, 
which means that we can take 7at+i := Hsl > 0. □ 

Now the assertion of part (ii) of the theorem is an immediate consequence 
of Lemmas 2 and 3. Theorem 1 is proved. □ 
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3 Local behaviour — a crude explanatory model 



Simulations of the f-process with S(v) = v a show avalanche-scale clustering 
leading to the formation of a permanent tight cluster when a > 1 and N is 
sufficiently large, and a weak variable clustering when a < 1. It would be 
of interest to obtain some insight into the causes of this phase-change type 
of phenomenon. In this subsection we will present a simplistic model at a 
physicist's level of rigour which explains why such a transition occurs at the 
threshold value a — 1. 

Let A C M be some small connected "test region" and Na(x) := card({i : Xi G 
We consider the evolution of ]Va(X„). Let 

y B {x) :^S(^pcM. 

Given X n = x, at the following step of the process the probability of a point 
being lost from A is 

sUxj' (18) 

while the probability of a fresh point entering A is n(A). As the total number 
iV of points in the process is typically much larger than N A (X. n ), it is natural to 
expect that, due to the effect of a law of large numbers, the relative fluctuations 
in 5?m (X n ) will be relatively small compared to those in S^aO^ti) (this is borne 
out by the results of simulations). So let's assume, for simplicity, that 5? = 
y M (y^ n )) is constant. 

We wish to find an approximation to (|18p as a function of N A = N A (X n ). 
Another simplifying assumption (also supported by simulations) is that the 
conditional (given Na) distribution of cell volumes for cells whose generators 
lie in A is the same, modulo scale, for different values of Na- That is, for the 
conditional (on Na) distribution function of the volume V of a randomly chosen 
cell with generator in A one has 

Pr(V <v\N A )=g(v/m A ), v > 0, 
for some g, where uia = ~Ei(V\Na)- Furthermore, we can take 

m A = N A - 1 E HC%:JNa)*{3\(A)/N a , (19) 

where /3 is a quantity dependent only on the geometry of A and the order 
of magnitude of N A , which reflects the fact that the union of the cells with 
generators in A overlaps A itself. A rough calculation shows that (3 = 1 + 
0(N^ 1 ). When A is a square or a circle and Na is not often less than 20, (3 
can be considered to be in the range (1.0, 1.5]. Appealing to the law of large 
numbers, we could write 

y A {^n) « NaE (S (A(C^" j )) I Na) = N A / S{v)dg{v/m A ), 

Jo 
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where / stands for the index of a "typical" cell with a generator X n j £ A. Now 
since S'(u) = v a , using (JT§j) the integral above becomes 

p\(M) p\(M)/m A 

/ v a dg(v/m A ) = m A / w a dg(w) 
Jo Jo 

POO 

« iVr(/3A(M)) Q / 
./o 

making the natural assumption that the last integral converges. Combining 
these approximations we get the following estimate for the probability that a 
point is removed from A in one step: 

POO 

^A^y- 1 « N^' a ((3X(M)) a y~ 1 / ^(ty) = KN\- a , 

Jo 

for some constant K. Thus, if AN A denotes the change in N A in one step of 
the v-process, we have 

E(AN A \N A ) w - KN\~ a . (20) 

As the right-hand side of (f2*Uj) is an increasing function of iV^ when a > 1, 
in this case we have a positive feedback condition for the mean of the number 
Na of points in our test region A. This means that the process is bound to 
quickly leave the "intermediate" range of states characterised by diffuse, roughly 
uniform spatial point distributions — an observation that is in agreement with 
the simulation data. Note also that, once the "destabilising mechanisms" have 
transformed the point distribution to a single (or a few) tight cluster (s), the 
assumptions on which the crude local model was based are no longer valid. 

On the other hand, when a < 1, the relation ()2()j) expresses either neutral 
(a = 1) or negative feedback. Hence one expects "local stability" from the 
process behaviour: small clusters of points will form and disappear without any 
"global" dramatic changes for the whole picture. 

In such cases we can expect N A to take values close to Nfi(A), so we can 
estimate K by the equation 

K(Nfj,(A)) 1 " a w fi(A), 
from which we conclude that 

K ss fi(A) a N a -\ a<l. 
This approximation is reasonably well-supported by simulations. 
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